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ABSTRACT 

We investigate the formation of bow shocks around exoplanets as a result of the interaction 
of the planet with the coronal material of the host star, focusing on physical causes that can 
lead to temporal variations in the shock characteristics. We recently suggested that WASP- 
12b may host a bow shock around its magnetosphere, similarly to the one observed around 
the Earth. For WASP 12b, the shock is detected in the near-UV transit light curve. Observa- 
tional follow-up suggests that the near-UV light curve presents temporal variations, which 
may indicate that the stand-off distance between the shock and the planet is varying. This 
implies that the size of the planet's magnetosphere is adjusting itself in response to variations 
in the surrounding ambient medium. We investigate possible causes of shock variations for 
the known eccentric (e > 0.3) transiting planets. We show that, because the distance from the 
star changes along the orbit of an eccentric planet, the shock characteristics are modulated 
by orbital phase. During phases where the planet lies inside (outside) the corotation radius of 
its host star, shock is formed ahead of (behind) the planetary motion. We predict time offsets 
between the beginnings of the near-UV and optical light curves that are, in general, less than 
the transit duration. Variations in shock characteristics caused in eccentric systems can only 
be probed if the shock is observed at different orbital phases, which is, in general, not the case 
for transit observations. However, non-thermal radio emission produced by the interaction of 
the star and planet should be modulated by orbital phase. We also quantify the response of 
the shock to variations in the coronal material itself due to, e.g., a non-axisymmetric stellar 
corona, planetary obliquity (which may allow the planet to move through different regions 
of the host star's corona), intrinsic variations of the stellar magnetic field (resulting in stellar 
wind changes, coronal mass ejections, magnetic cycles). Such variations do not depend on the 
system eccentricity. We conclude that, for systems where a shock is detectable through transit 
light curve observations, shock variations should be a common occurrence. 

Key words: planet-star interactions — planets and satellites: individual (HD 80606b, HD 
17156b, CoRoT-lOb, HAT-P-2b, HAT-P-17b, WASP-8b) — planets and satellites: magnetic 
fields — stars: coronae — stars: winds, outflows 



1 INTRODUCTION 

More than 100 exoplanets have now been detected through optical 
transit observations. Transit light curves can provide a wealth of 
information about both the object in transit and its host star, such 
as the planet radius R,,, orbital pe riod Ppih, orbital semi-major axi s 
a, stellar mass M, and radius R, jSeager & Mallen-Omelasll2003h . 
Observations of light curves in wavelengths other than in the op- 
tical should also provide more information about the system. This 
is the case of the gas giant planet WASP-12b, whose near-UV light 
curve differs from t he optical one by presenting an earlier beginning 
jpossati et alj2O10l) . This early ingress has been interpreted as due 
to the presence of material around the planet, although several ex- 
planations for the origin of this material have been suggested in the 
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literature jLai et al.ll2010l: I Vidotto. Jardine & HellindlloTol . here- 
after Paper 1). In Paper ll we suggested that this material probes 
the existence of a bow shock around the planet, which is formed 
when the relative motion between the planet and the stellar coronal 
material is supersonic. A bow shock is observed around the Earth's 
magnetosphere as a result of the interaction of our planet with the 
solar coronal wind. According to our model, the detection of the 
shock through transit observations can constrain the planetary mag- 
netic field, something that is as yet unobserved in exoplanets. 

In Ividotto, Jardine & Hellin hereafter Paper 2), we 

applied the shock model initially developed for WASP-12b to the 
transiting systems known at that point, determining which planets 
are prone to develop shocks. We predicted that a significant number 
of transiting systems (36 out of 92 planets) might have a detectable 
shock, implying that bow shocks might indeed be a common fea- 
ture surrounding transiting planets. These shocks can cause light 
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curve asymmetries, that once observed can allow us to determine 
the planetary magnetic field intensity . 

In our previous works /P aper lUPaper 2l) . we considered these 
shocks to be static. However, preliminary results from a second set 
of near-UV observations from Haswell et al. (in prep.) show that 
the early ingress observed in the near-UV light curve of WASP- 12b 
presents temporal variations. They found that the early ingress was 
more delayed in the observations of Apr/2010 than in the first set 
of observations from Sep/2009. This has motivated us to investigate 
which effects could cause temporal variation in the near-UV light 
curve of a transiting planet. 

Such temporal variation in the light curve of WASP-12b indi- 
cates that the stand-off distance between the shock and the planet is 
varying. According to our model, the shock limits the extension of 
the magnetosphere of the planet, so that variation in the stand-off 
distance implies that the size of the planet's magnetosphere is ad- 
justing itself in response to variations in the surrounding ambient 
medium. 

In this paper we consider possible causes of time-dependence 
in the location and density of the shocked material around the 
planet. This will occur if the planet experiences changes in its lo- 
cal environment (such as the density, magnetic field strength and 
stellar wind speed) as it moves in its orbit. Three causes might be: 

(i) The stellar magnetic field is in a steady state and is az- 
imuthally symmetric, but the planet is in an eccentric orbit. There- 
fore, during the orbit, the distance from the planet to the host star 
varies. 

(ii) The stellar magnetic field is in a steady state but is not az- 
imuthally symmetric. During its orbit the planet therefore experi- 
ences different local conditions due to these azimuthal magnetic 
field variations. 

(iii) The stellar magnetic field is intrinsically time-dependent - 
either due to turbulent fluctuations in the stellar wind, flaring, coro- 
nal mass ejections, or longer-timescale variation due to a magnetic 
cycle. 

Our aim is to investigate how each of these processes may 
contribute to temporal variation on the planetary bow shock and 
magnetosphere. In this paper, we cons ider only systems with high 
eccentricities, which were neglected in lPaper 2| 

This paper is organised as follows: Section |2] discusses the 
effect on the shock configuration around a planet in an eccen- 
tric orbit. Section |3] investigates the three suggested causes of 
time-dependence in light curve asymmetries. Because the distance 
from the star changes along the orbit of an eccentric planet, the 
shock characteristics are modulated by orbital phase. Similarly, 
non-thermal planetary radio emission should be modulated by or- 
bital phase. We estimate the amplitude of the modulation on the 
frequency of radio emission in Section |3] A summary and conclu- 
sions are present in Section|4] 



2 BOW SHOCKS AROUND ECCENTRIC PLANETS 

In order for a bow shock to form, the condition that needs to be 
satisfied is that the planetary velocity through the stellar coronal 
material has to be supersonic. We investigate two different scenar- 
ios for the stellar coronal material, which can be (1) static, as in a 
corona that is magnetically confined to corotate with the host star 
or (2) dynamic, as a result of the presence of a stellar wind. Shocks 
can be developed in both scenarios, but the orientation of the bow 
shock (given by the angle between the normal of the shock and the 
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Figure 1. Sketch of shock formation. For positions in the orbit where the 
planet-star distance nr is nr < rco, an ahead-shock is foimed. For -m > r^o, a 
behind-shock exists. 



azimuthal vector - see Fig. 1 in IPaper ih varies depending on the 
scenario adopted. In scenario (1), the angle is either 6 = 0°, when 
the shock forms ahead of the planetary orbital motion, or 9 = 180°, 
when the shock forms behind the planetary orbital motion. In sce- 
nario (2), this angle lies between 9 = 0° and 90°, in general. 

The condition that defines whether a shock formed in the mag- 
netically confined case is an "ahead-shock" (9 = 0°) or a "behind- 
shock" (6 = 180°) depends on whether the coronal material rotates 
faster than the orbital velocity of the planet. If a prograde planet 
orbits the star beyond its Keplerian corotation radius, the coronal 
plasma lags behind the planetary motion, and the shocked material 
accumulates behind the planetary orbital path, trailing the planet. 
This gives rise to a behind-shock. CoRoT-Ilb is a candidate to 
host a shock in such a configuration. In the opposite case when the 
planet orbits the star inside its corotation radius, the shock forms 
ahead of the planetary motion. 

We note that certain eccentric systems can cross the corotation 
radius twice during their orbits (Fig.[T|. The corotation radius r^o is 
given by 



GM. 



R,[v s'm(i)/ sin(0] 



1/3 



(1) 



During the orbital phases where the distance between the planet 
and the star (ru) lies inside rco, the shock forms ahead of the plan- 
etary motion. For the remaining phases of the orbit (m > r^o), the 
shock trails the planet. At the two points in the orbit where ru = r^a, 
the shock weakens considerably and eventually may disappear, be- 
cause the net azimuthal velocity between the coronal material and 
the motion of the planet is small. 

Table [T] presents the planetary and stellar parameters for 
the most eccentric transiting planets (e > 0.3) known to date. 
The data used here are taken from the exoplanets encyclopedia 
(http://exoplanet.eu/catalog.php I, except for the sky-projected stel- 
lar rotati on veloci ty v sin(i), which was taken from the compilation 
made by S chlaufmanI ( 120 1 01) for most of the systems . For CoRoT- 
10 and HA T-P-I7, values of v si n(i) were taken from ^Bonomo et al] 
( I2OIOI) and lHoward et"ai] ( l20IOl) . respectively. Values of impact pa- 
rameter b and transit duration t^^ were compiled fro m the litera- 
turj] (LM outou et al. 2009; Winn et al. 2009; Bonom o et alj|2010l : 
[Howard et alj20IOl : iPal et al.ll20ia ; lOueloz et alj|201(]l) ^ 

For an elliptical orbit, the distance ru from the planet to the 
star (located at one of the focus of the ellipse) is given by 



For the systems b was not available, we calculated it as b 
(a/R,)cos!(l -e^)/(l + e sin w), where ai is the argument at periastron. 
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Table 1. Planetary and stellar parameters for the most eccentric transiting planets (e > 0.3) known to date. The columns are: (1) planet name, (2) planet mass, 
(3) planet radius, (4) planet orbital period, (5) eccentricity of the orbit, (6) semi-major axis, (7) optical transit duration, (8) impact parameter, distance to the 
star when the planet is (7) at periastron and (8) at apastron, (9) stellar mass, (10) stellar radius, (11) sky-projected rotational velocity, (12) stellar corotation 
radius, and (13) the ratio of the size of planet magnetosphere at apastron and at periastron. 
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where e is the eccentricity and is the angle formed between the 
connecting planet-star line at nearest approach (periastron) and the 
position of the planet. Therefore, at periastron, = 0°, while at 
apastron = 180°. The radial and tangential velocities of the 
planet (assuming M, » Mp) are, respectively, 
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[1 -l-ecos(0)]. 
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The condition for a shock to be formed is that the relative ve- 
locity between the material surrounding the planet and the planet 
itself |Au| must exceed the sound speed c, = (kgT/my^^, where 
ks is the Boltzmann constant, m = /imp is the mean particle mass, 
with fi = 0.66 and m,, the proton mass. Therefore, there exists a 
maximum temperature below which shock formation is allowed 



lAul-ni 



(5) 



Here, we make a distinction between two possible scenarios: 
(1) a corona that is magnetically confined to corotate with the host 
star and (2) a case where the coronal material expands in the form 
of a stellar wind. In the magnetically confined scenario, the coronal 
material has only rotation (azimuthal) velocity given by 



VSin(i) T77 

sin(;) R, ' 



(6) 



where the angular velocity of the star Q., = [vsin(!)/ sin(i)]/R,. In 
the wind scenario, the material surrounding the planet has both a 
rotation velocity (conservation of angular momentum of particles 
in the wind) 



£l,Rl V sin(0 R, 



sin(i) ru 



(7) 



and a radial velocity Mr,,,, due to the stellar wind. By as suming 
that t he host star has an isothermal, thermally-driven wind jParkeJ 
Il958 f). we calculate the wind velocity profile from the integration 
of the differential equation 



or 



' dr 



GM, 



(8) 



where n„ is the wind density, r is the radial distance from the star, 
p,„ = rin.kgT is the wind thermal pressure. 

Therefore, the relative velocity |Au| between the material sur- 
rounding the planet and the planet itself for the confined case is 

|AU| = lUcor - Vp,| = [v;_p, + (li^,eor - %?/)']"' (9) 

and for the wind case is 



|AU| = |U„ - Vp/I = [(m,.„ - Vr,pl) + - U^,plY\ 



2.1 HAT-P-2b 



211/2 



(10) 



HAT-P-2b is a good illustrative example to compare both scenarios 
for the suiTounding stellar material, since its orbit may lie almost 
entirely within the stellar magnetosphere and, in this case, the hy- 
pothesis adopted in the confined scenario case should be met during 
almost its entire orbit. From equations (|5), (|2) and ( 1101 . we show 
in Fig. [2] the maximum temperature required for shock formation 
for HAT-P-2b for both the confined corona scenario (black) and the 
wind scenario (red), where we adopted a coronal temperature of 
r = 2 X 10* K in the evaluation of u,.„. (equation[8). It is interesting 
to note that at orbital phases near ^ 101° and 259°, the 
velocity difference of the confined corona case approaches its min- 
imum and a shock will only be formed if the temperature is very 
small (r < Tniiix - 5 X 10^ K). Because such a small temperature is 
not expected in the corona of a F8-star, it is probable that a shock 
will not be formed near these orbital phases. Furthermore, for the 
confined corona scenario, if the coronal temperature of HAT-P-2b 
is r = 2 X 10* K, a shock will not be formed for orbital phases 
where r„,„^ > T = 2 x 10* K, i.e., between ^ 158° and 212° 
(around apastron). 

We note that even if a shock is formed around the planet, 
it may not be detected if the compressed local plasma does not 
achieve a density high enough to cause an observable level of op- 
tical depth. For the corona case, the local density of the corona is 
given by 



: riQ exp 



GM,/R, 
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^orb 
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where no is the density at the base of the corona. For simplicity, 
we assume that no = 10* cm"^, similar to that of the solar corona 
(Withbroe 1988). For the wind case, the density n„, is given by con- 
servation of mass of the wind n„,w,,„,r^ = constant, with Ur,„. given 
by equation l[8j. Fig. |3] shows how the density varies according 
to the planetary orbital phase (given by 0) for HAT-P-2b assum- 
ing two difi"erent temperatures for the material that surrounds the 
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Figure 2. Maximum temperature required for sliock formation for HAT-P- 
2b for botli the confined corona scenario (black) and the wind scenario (red). 
At = 0, 360°, the planet is at periastron and at = 180° at apastron. At 
orbital phases corresponding to =: 101° and ^ 259°, the net velocity 
between the planet and the coronal material is very small, so the temperature 
required for shock formation in the confined corona scenario should also 
be very small, implying that a shock will not be formed near these orbital 
phases. Considering T = 2x 10* K for the wind scenario. 

planet. Because HAT-P-2 has a relatively high v sin(;), the density 
at apastron is larger than at periastron in the confined case (the sec- 
ond term in the RHS of equation dill ) dominates over the first term). 
We note, however, that these 'bumps' may not be physical, since the 
corona may not still be corotating with the star at large distances. A 
way to estimate where this assumption breaks down is to evaluate 
the plasma-yS, the ratio between thermal and magnetic energy densi- 
ties: whenever the thermal pressure takes over the plasma pressure, 
the hyp othesis may become invalidated. 

In|Paper 2, we suggested that a local density surrounding the 
planet as low as 10* cm"^ can still provide detection. Therefore, 
according to Fig. [3] detection of shocks around HAT-P-2b implies 
that the coronal temperature should be T > 2 x 10* K. Note, how- 
ever, that because a shock does not exist during orbital phases 
158° < < 212° (when T > r„,„, assuming T = 2 x 10* K), 
shock formation and detection around a planet can constrain both 
its surrounding density and temperature. 

We define the angle that the shock normal makes to the relative 
azimuthal velocity of the planet as 

e = arctan( ^ \ (12) 

for the confined corona scenario and for the wind scenario as 

l\v,.pl-Ur,M,\\ 

e = arctan . (13) 

Due to the ©-dependence of the involved velocities, the angle of 
the shock 6 also depends on the orbital phase. Fig. |4] shows the 
shock angle d for both the corona (black) and wind (red) scenarios 
for HAT-P-2b assuming T = 2 x 10* K. For eccentric systems in 
the confined scenario, we note that a purely ahead-shock (6 = 0°) or 
behind-shock (9 = 180°) only occurs when the radial component of 
the orbital velocity is null (v^j,/ = in equation [T2}, which happens 
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en 

Figure 3. Density of the ambient medium through where the planet HAT-P- 
2b orbits as a function of planetary orbital phase (given by 0). Two dift'erent 
temperatures are assumed: T = 1 X 10* K (dashed lines) and 7 = 2 X 10* K 
(solid lines), for both the confined corona (black lines) and the wind (red 
lines) scenarios. 




Figure 4. The angle 6 that the shock normal makes to the relative azimuthal 
velocity of the planet as a function of orbital phase given by 0, for HAT- 
P-2b. A purely ahead-shock (0 = 0°) occurs at periastron (0 = 0,360°). 
A purely behind-shock {9 = 180°) occurs at apastron (0 = 180°). A day- 
side shock (8 = 90°) occurs at = 98, 262° (positions marked with filled 
circles). 



at periastron (0 = 0,360°) and apastron (0 = 180°), respectively. 
A purely day-side shock (9 = 90°) happens when «^_„,r - v^ ,,; = 0, 
which occurs at orbital phases = 98, 262°, indicated by filled 
circles in Fig.|4] For the wind case, because m^ ,,,, < v^ pi in general, 
9 < 90°, implying that the shock never forms behind the planetary 
motion. For hig her wind temperatures, we expect the shock angle 
to approach 90° jPaper ih . 
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The extent of the planetary magnetosphere can be deter 
mined by static pressure balancq^ 



, 2 B{mf 
innAu -\ 1- p 



8;r 



-Pp- 



(14) 



Both n and p are the local density and pressure of the material sur- 
rounding the planet, which can either refer to the coronal scenario 
("cori Pcoi = ricoiksT) or to the wind scenario (n„„ p„,). B(m) is the 
stellar magnetic field intensity at the site of the interaction with the 
planet. In the confined case, we assume that the stellar magnetic 
field is dipolar so that B{m) = B,(R,/m)^. In the case where the 
corona escapes in the form of a wind, we neglect the effect of the 
stellar magnetic field. We also consider that the planet magnetic 
field geometry can be approximated as that of a dipolar field, so 
Bpir/n) = BpiRp/r/n f, where Bp and Bpirm) are the magnetic field 
strengths at the surface of the planet and at the nose of the planet's 
magnetosphere, respectively. We neglect the planet thermal pres- 
sure Pp. 

If the dominant terms in equation ( I14t are the magnetic terms, 
we showed that, for the confined corona scenario dPaper ll) , equa- 
tion l ll4t reduces to 



rM_ 

Rn 



<B 



"3 ^ 

S. I R,' 



(15) 



In this case, the ratio between at apastron and at periastron is 
»"M(?3Tap) 1 -I- e 



'■M(ra-pc) 1 



(16) 



where we used equation (|2). For HAT-P-2b, rMi:^^,)! rmi'^j^pi) ~ 
3.1. Table[T] column 13 presents this value for the remaining plan- 
ets. 

In our previous works dPaper ll: IPaper 2b . for the confined 
corona case with a planet in a circular orbit, we assumed that the 
observed stand-off distance of the shock from the planetary sur- 
face is approximately the extent of the planetary magnetosphere 
tm- The stand-off distance can be observed through the time dif- 
ference of th e ingress of two transit light curves taken at different 
wavelengths jPaper ih . Therefore, it is straightforward to derive the 
size of the magnetosphere of the planet and the ratio Bp/B, for cir- 
cular orbits in a confined corona. If the magnetic field of the star 
is known, then the planetary magnetic field is easily derived. For 
the other cases (i.e., in the wind scenario and/or for eccentric sys- 
tems), the shock is not always formed purely ahead of the planet, 
as shown in Fig.|4] Therefore, from observations, one can derived 
the stand-off distance of the shock that is projected on the plane of 
the sky. 

For the calculations presented hereafter, we consider all the 
terms in equation il4i , except pp. Detailed calculation of pres- 
sure balance as a means to determine the size of a planet's mag- 
netosphere ha s been used by the Solar Sy stem community for sev- 
eral decades jChapman & Ferrarolll930h . In these calculations, a 
correcting factor is often assumed in the value of Bp(rM) (equa- 
tion [14) to account for the compression of the magnetic field at 
the nose of the magnetosphere. In the case of a dipolar planetary 
magnetic field, the presence of magnetopause currents gives rise 
to a magnetic field with the same i ntensity as the magnetic field 
just inside the magnetosphere (e.g., ICraven sl ll997UMever-Vemetl 



^ We are correcting a typographical error in equation (9) of ! Paper ll and 
equation (5) of lPaper 2. This has no effect on the calculations done in these 
papers. 
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Figure 5. The magnetospheric size rM of HAT-P-2b as a function of plane- 
tary orbital phase given by assuming a coronal temperature T = 2x 10* K 



and B„ 



14 G. Black curves refer to the confined corona scenario with dif- 



ferent 6, and the wind scenario is represented by the red solid line. 



doubling, therefore, the value of Bpiru)- To allow for dif- 
ferent magnetopause geometries, it is also considered that the mag- 
netic fiel d at the nose of the m a gnetopause increases by a factor 
/ (e.g., iGrieBmeieretaLl liioi iKhodachenko et aD l2007h . Alto- 
gether, these modifications result that BpirM) — > IfBpir^) in equa- 
tion dl4t . For the Earth's magnetospher e, the compression factor i s 
empirically determined as 2/ 2.44 jKivelson & Russelil 19951 ). 
Given the other uncertainties involved in the calculation of for 
exoplanets, we do not include the compression factor 2/ in this 
work. If it were considered, the magnetospheric radii derived in 
this paper would have increased by a factor of (2/)'^'. 

Fig. [5] shows the magnetospheric radius of HAT-P-2b derived 
using equation | |14| | for both the wind case (red) and the confined 
corona case (black) for Bp = 14 G (similar to that of Jupiter) and 
S. = 1 G (solid), 10 G (dashed) and 100 G (dot-dashed). For the 
confined case, we note that the more intense B, is, the more com- 
pact is the magnetosphere of HAT-P-2b. The magnetospheric sizes 
of the Earth, Jupiter, Saturn, Uranus and Neptune are 11, 5 0-100, 
16 - 22, 18, and 23 - 26 planetary radii, respectively dBagenall 
I1992I) . These values are larger than the values derived for HAT- 
P-2b (Fig.jSj- The compact size of the magnetosphere of HAT-P-2b 
is because the planet is located much closer to the star where the 
stellar wind pressures are still significant. 

The time difference 6t between the beginning of the optical 
and near-UV transits led us to sugge st the ex istence of a bow shock 
around the hot-Jupiter WASP-12b dPaper lb . This time difference 
can be infeiTed by geometrical considerations of the transiting sys- 
tems. Consider the sketches presented in Fig. |6] and d\}\i are, 
respectively, the sky-projected distances that the planet (optical) 
and the system planet-l-magnetosphere (near-UV) travel from the 
beginning of the transit until the middle of the optical transit 



and 



duv = {Rl - + rM, 



(17) 



(18) 



where b is the impact parameter derived from transit observations. 
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Figure 6. Sketches of the hght curves obtained through observations in the 
a) optical and b) near-UV, where the bow shock sun'ounding the planet's 
magnetosphere is also able to absorb stellar radiation. 



The optical transit duration fji is the time the planet takes to travel 
2dop (from beginning to end of the optical transit). Therefore, by a 
linear extrapolation we derive 

Adopting b = 0.395 R, and ta, = 0.1787 d jPal et alJIlOlOl) . Fig.|7] 
shows the time-offset 6t for HAT-P-2b for the confined corona 
(black) and wind (red) scenarios. We adopt Bp = 14 G and B, = 1 G 
(solid), 10 G (dashed) and 100 G (dot-dashed). We again note the 
modulation of St with planetary orbital phase. For B, > 10 G, 6t 
becomes comparable do the time offset observed for WASP- 12b 
(6t ^ 0.5 h). 



2.2 Remaining Eccentric Systems 

For a shock to exist and be observable, the following conditions 
must be met: 1) the temperature of the ambient medium surround- 
ing the planet has to be T < r^ax (equation O; 2) the local density 
of this ambient has to be larger than a minimum threshold to pro- 
vide enough detectable absorption (n> 10** cm"^); 3) the magneto- 
sphere of the planet has to be large enough to provide a time offset 
that is larger than the instrumental uncertainties in time. 

In order to evaluate the time differences St between the begin- 
nings of the optical and near-UV transits for the remaining systems 
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Figure 7. Expected time differences St between the beginnings of the optical 
and near-UV transits calculated for HAT-P-2b. Curve labels are as in Fig.|5] 
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Figure 8. Time-delay St between optical and near-UV transit beginnings 
scaled to the transit duration ?(], for the eccentric systems studied in this pa- 
per (Table[T]. Planets are shown in order of increasing eccentricity: WASP- 
8b, HAT-P-17b, HAT-P-2b, CoRoT-lOb, HD 17156b, HD80606b. 



presented in Table [T] we performed a similar analysis to that pre- 
sented in Section im We assume that all host stars analysed in this 
paper have coronal winds with T = 2 x 10' K and base densities 
no = 10** cm"^ and planets have B,, = 14 G. We evaluated St for 
both the confined corona and wind cases. Fig. [8] shows the min- 
imum and maximum expected time offset St scaled to the transit 
duration fjr for a stellar magnetic field of fi, = 1 G. The confined 
corona scenario (black) provides smaller time offsets than the ones 
derived in the wind scenario (red). In Fig.|8] we assume that shocks 
can be formed at all points in the planetary orbit, although this may 
not always be the case (see previous paragraph). 

Near-UV light cu rve asymmetry has been detected in the gi- 
ant planet WASP- 12b l lFossati et al.ll201(il) . So far, we do not know 
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the precise location of the material that is causing the absorption 
observed in Mg II lines. It may be that absorption is concentrated 
in a certain location of the planetary atmosphere (e.g., only around 
the shocked material) or that it extends through a larger area of the 
planetary atmosphere. Observations of near-UV light curves from 
different systems can help constrain the location of absorption. 

Consider the case where the material causing absorption in 
Mg II lines is restricted around the small region of the shocked ma- 
terial, which is located at a distance ~ vm from the center of the 
planet. If this distance is larger than the path length of the tran- 
sit (2[/fJ - fc^]''^), two independent transits occurs (Fig. |9). The 
first one is the transit of the shocked material itself, which ends 
before the planetary transit begins. After the first transit, the light 
curve returns to the continuum level until the planetary transit it- 
self takes place. We note that if one just observes immediately be- 
fore the planetary transit, the transit of the shocked material may 
be missed. Now, consider that the material causing absorption in 
Mg II lines fills the magnetosphere of the planet, still assuming that 
tm > 2{Rl-b^yi^ . Instead of causing a double transit, the transit of 
the shocked material will blend with the transit of the planet, and 
only one transit occurs. 

The condition > 2{R\ - hP-)^^^ can be written in terms of 
the transit offset and duration as Stlt^, > I (Stlti^ = I: dashed line 
in Fig. [8}. Among the eccentric systems analysed in this paper with 
the assumed values of To = 2 x lO*" K, Bp = 14 G, S. = 1 G 
and iiQ = 10** cm"', we note that there are only two systems suit- 
able to probe the extension of the absorbing material: HD 80606b 
and CoRoTlOb. These systems could present a double transit if ab- 
sorption is restricted to a small area surrounding the shock or a 
single transit if absorption occurs at an extended area. We note that 
the largest values of time offsets in Fig.[8]happen during apastron, 
when the size of the magnetosphere is the largest one (wind case). 
However, during apastron the material surrounding the planet has 
its smallest density. Therefore, it is more difficult to observed shock 
at apastron, where double transits are more probable to happen. 

The case where the magnetosphere of the planet is contained 
inside the projected distance (i.e., r/n < (Rl-b^Y'^ or, equivalently, 
dt < f(ii) always results in a single transit. 



3 DISCUSSION 

Bow shocks around planets may not be steady. If, for instance, a 
planet has a thermal mass loss that intrinsically varies, the bow 
shock formed at the edge of the planetary magnetosphere could 
expand whenever there is any increase in the planetary mass-loss 
rat e (for mass-loss v ariations due to the evolution of the host star, 
see lMurrav-Clav eTal. 2009) . Also, if the dynamo responsible for 
generating the planetary magnetic field changes its efficiency, the 
size of the planet's magnetosphere should change accordingly. 

Another circumstance that may cause the position and den- 
sity of the shocked material around the planet to vary is related to 
variations not in the planet itself, but in the ambient medium that 
surrounds it. The last alternative is the focus of our paper. Here, we 
discuss three possible causes that could lead to such variations. 



3.1 Case 1: Eccentric Planet 

The first case we discuss is the case where the confined corona or 
coronal wind is in a steady-state and is axisymmetric. We also as- 
sume that the orbital axis of the planet is aligned with the stellar 



axis of symmetry. If the planet is in a circular orbit, during its revo- 
lution around the star, the stellar material encountered by the planet 
does not change its properties. Therefore, in the reference frame of 
the planet, the shock is steady. However, an eccentric planet probes 
different plasma conditions as its distance to the host star varies 
through its orbit. This case was detailed in Section |2] where we 
showed that the characteristics of the bow shock formed around the 
planet is modulated by the orbital phase. We also showed that at 
certain phases, a shock may not be formed (depending on the tem- 
perature of the coronal material) and that even if a shock is formed, 
it may not have enough density to be detected. 

This variation in shock characteristics can only be probed if 
the shock is observed at different orbital phases, which is not the 
case for transit observations that are, in general, conducted at only 
one orbital phase. However, any radio emission produced by the 
interaction of the star and planet should be modulated by orbital 
phase. The amplitude of such modulation is estimated next. 



3.1.1 Planetary Radio Emission 

The four giant planets of the Solar System and the Earth emit at ra- 
dio wavelengths. Because there are correlations between this emis- 
sion and the local characteristics of the solar wind, planetary ra- 
dio emission is associated to the interaction of the solar wind with 
the planets. In analogy to the Solar System, it has been suggested 
that exoplanets may also produce radio emission when they i nteract 
with the winds of their host stars jjardine & CameronlllOOsl) . 

In Jupiter, decametric emission is thought to arise in a ring 
surrounding the auroral region wherein the planetary magnetic field 
lines are open. The aperture of the auroral ring can be related to the 
size of the planet's magnetosphere rM in the following way. We 
consider the planet to have a dipolar magnetic field aligned with 
the planetary orbital spin axis, such that 



B„ 



and 



B„ IR 



B,,.a = — I — I sintt. 



(20) 



(21) 



where R is the radial distance from the planet. Bp is the planetary 
magnetic field evaluated at the pole and a is the angle that mea- 
sures the colatitude (Fig.lIO). The paths of individual field lines are 
described by A = sin" alR where the value of A labels each field 
line (i.e., it is constant for each line). If the largest closed field line 
connects the surface at a colatitude q-q (where A = sin' ag/Rp) to 
an equatorial radius R = vm (where A = 1 /r^), then 



1 sin" cfo 



(22) 



The fractional area of the planetary surface that has open magnetic 
field lines is then (I - cos uq), where we considered both the north 
and south auroral caps of the planet. Thus, for example, if vm = 
5 Rp, then q-q - 27° and open field covers ~ 11% of the surface. 
If = 2 Rp, then ao - 45° and open field covers ~ 29% of the 
surface. If the planetary radio emission comes from the region that 
separates closed and open field lines (auroral ring), such emission 
would happen in a cone with aperture given by 



ao = arcsin [(Rp/r„) '-]. 

The planetary magnetic field at this colatitude ao is 



(23) 
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Figure 9. Sltetch of the double transit caused by a very extended planetary magnetosphere. If the material absorbing near-UV stellar radiation is concentrated 
in a small region around the bow shock (see shaded area), it ingresses in transit before the planetary transit. This causes the first dip in the light curve. After its 
transit, the light curve returns to its continuum level until the planet itself ingresses in the transit (second dip). If the absorbing region fills the magnetosphere 
of the planet, the fi'ansit of the material and the planet are blended together and only one transit is observed. 




Figure 10. Sketch of the planet's magnetic field lines. The closed-field line 
region extends out to tm ■ The area where inagnetic field lines are open and 
exposed to the stellar wind (auroral region) extends out to planetary colat- 
itudes ffo, which is determined by the magnetosphere 's size ru (equation 
123) . Radio emission is believed to aiise from a ring sunounding the auroral 
region. 



B(ao) = ^(1 + 3 cos^ ao)"', (24) 
which corresponds to an electron cyclotron frequency 

/, = 2.8|^^J MHz. (25) 

Despite several searches, radio emission from exoplanets has 
not yet been detected. One possible reason for that may be that 
emission occurs at a frequency much different from the frequency 
range of the radio detectors used in the observations. The frequency 
of the radio emission is believed to be a fraction of /,.. As the plan- 
etary magnetospheric size is modulated by its orbital phase in an 
eccentric system (Section|2}, the latitude of the last closed field line 
and, consequently, the fraction of the area of the planet that will be 
exposed to the stellar wind are also modulated by the orbital phase. 
Therefore, we expect that will present a similar modulation and 
the frequency where radio emission occurs should differ at difFer- 
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Figure 11. Cyclotron frequency fc of the electrons believed to be respon- 
sible for generating planetary radio emission. The bandwidth where radio 
emission should occurs believed to be a fraction of . Planets are shown in 
order of increasing eccentricity: WASP-8b, HAT-P-I7b, HAT-P-2b, CoRoT- 
10b, HDI7I56b, HD80606b. 



ent parts of the orbit. In Fig.[TT] we present the minimum and max- 
imum cyclotron frequencies for each of the systems presented in 
Tablefflfor the case S,, = 14 G and B. = 1 G. We note from Fig.fTTI 
that fc is larger for the wind scenario than for the confined case. 
This is because the size of the planetary magnetosphere is smaller 
for the latter one. We also note that the modulation in fc does not 
vary significantly with planetary orbital phase, reaching an ampli- 
tude of only a few MHz and a frequency centred at around 37 MHz. 
Although small, this modulation could be spectrally resolved with 
radio telescopes, such as LOFAR. 

Because radio emission is related to the power of the impact- 
ing coronal material on the planet, it is preferable to look for radio 
emission during the closest approach of the planet. Radio emission 
from HD80606b d uring the periastron was recently investigated by 
iLazio et alj i2010l) . although the detection was unsuccessful. Be- 
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cause the star-planet distance during periastron for the HAT-P-2b is 
about half the value of the same distance for HD80606b, HAT-P-2b 
should be a potential system to look for radio emission. 

3.2 Case 2: Azimuthal variations in the coronal material 

The second case we consider is the case where the coronal material 
is also in a steady state in the stellar reference frame, but it is non- 
axisymmetric. Because the stellar rotation is not in general locked 
to the planetary orbital period (except for the case the planet is in 
the corotation radius of the star), as the planet revolves around the 
star, it encounters material with different characteristics, causing 
shock variations. By performing several observations of transit, one 
can infer the impact on the planetary magnetosphere and bow shock 
of azimuthal variations in the stellar coronal material. 

An example of a non-axisymmetric coronal material is the 
case of an oblique stellar magnetosphere, where th e stellar spin 
axis is not parallel to the magnetic moment vector dVidotto et al.l 
lipid ). During its orbit, the planet may encounter a region of con- 
fined coronal material (e.g., when it is inside the closed magnetic 
loops of the host star), but at other orbital phases, it may be in the 
wind of the host star (where the magnetic field lines of the star 
are open and the stellar wind escapes). As the characteristics of 
the closed corona are different to the open-wind region as demon- 
strated in Section |2] (Figs. |2]-|7j, the magnetospheric size of the 
planet varies with orbital phase. So if, for instance, this is the case 
of HAT-P-2b, the time offset between optical and near-UV transits 
can span through a larger interval (e.g., 0.56 to 1. 1 h, for the condi- 
tions assumed in Fig.|7]and 6, = 1 G). Similarly, the variation in 
can be more significant than the variation we estimated individually 
for the wind scenario or the confined scenario. 

Shock variations as caused by azimuthal variations in the coro- 
nal material can be found in both circularised and eccentric orbits. 
We note that WASP- 12b (which is currently believed to be in a 
circular orbit) will also experience variations in the shock if the 
planet's orbit take s it thro ugh both regions of confined plasma and 
also wind plasma (Paper 1). In this case, azimuthal variations in the 
coronal material could explain the most recent results from Haswell 
et al. (in prep.), which show that the time-offset between the begin- 
nings of optical and near-UV transits varied from observations done 
in Sep/2009 and Apr/2010. 

We also note that, a similar effect could be observed in the 
case that the coronal material is axisymmetric (e.g., in an aligned 
stellar dipole), but the orbital plane and the stellar equator are not 
coplanar (non-null planetary obliquity). 

3.3 Case 3: Time-dependent Stellar Magnetic Field 

The last case we investigate in this paper is when the stellar mag- 
netic field is intrinsically time-dependent, which generates varia- 
tions in the stellar wind. Stellar magnetic field variations can oc- 
cur on different timescales, such as short timescales due to flares 
or coronal mass ejections (CMEs) and large ones as due to stellar 
magnetic cycles. 

When a CME hits a planet, it increases the density and velocity 
of the surrounding material, which may cause significant compres- 
sion in the size of the planet's magnetosphere. In this case, the den- 
sity contrast increases, which is favourable for shock detection. As 
the average duration of a solar CME at a distance 6-10 is about 
8 h jKhodachenko et alj|2007l) . the density increase effect caused 
by a CME is transient and may be easily missed, except for the 



case of planets orbiting young, magnetically active stars for which 
the magnitude and frequency of flares and coronal mass ejections 
may be much higher than for the present-day Sun. 

Although in the case of solar-like stars we may expect a mag- 
netic cycle on a timescale of decades (compared to 1 1 yr for the 
Sun), the rec ent detection of a 2-year cycle on the planet-hosting 
star T Bootis jpares et al.l2009t) may indicate that shorter cycles are 
possible. Therefore, variations in the magnetosphere of the planet 
due to polarity reversals of the stellar magnetic field should occur 
on timescales of magnetic cycles of the host star. 



4 SUMMARY AND CONCLUSIONS 

In this paper, we investigated time-dependent effects in the asym- 
metry of planetary transit light curves. As suggested by us, the 
presence of a planetary bow shock, formed from the interaction 
with the stellar coronal material, may be able to absorb stellar near- 
UV photons, causing the transit light curve in the near-UV to differ 
from the optical one. Our investigation of time-dependent effects in 
shock formation is motivated by preliminary results from a second 
set of near-UV observations of WASP- 12b obtained in Apr/2010 
(Haswell et al. in prep.). As compared to the first set from Sep/2009, 
it was found that the early ingress observed in Apr/2010 started be- 
fore the early ingress observed in Sep/2009. A possible explanation 
for this may be that the size of the planet's magnetosphere changed, 
therefore changing the location of the bow shock around the planet. 

We investigated three different scenarios that could cause the 
size of the planetary magnetosphere to change. In the first scenario, 
we considered the case of an axisymmetric stationary stellar corona 
and an eccentric planet whose orbital spin is parallel to the stellar 
axis of symmetry. Because of the planet's eccentricity, the coro- 
nal material interacting with the planet has different characteristics 
depending on the planetary orbital phase. Therefore, variations in 
the size of the planet's magnetosphere are modulated by the orbital 
phase. These variations can not be detected by observations con- 
ducted at one single orbital phase, which is, in general, the case for 
transit observations. 

However, observations that probe the size of the planet's mag- 
netosphere and that can be conducted at any orbital phase should 
be modulated by the orbital phase. This is the case for observations 
searching for non-thermal planetary radio emission due to the inter- 
action of the planet with the stellar material. This interaction limits 
the size of the planetary magnetosphere, which then constrains the 
size of the planet's area that is open to the stellar wind. We showed 
that the frequency of radio emission coming from this region will 
be modulated by the orbital phase, although its amplitude is at most 
a few MHz centred at ~ 37 MHz for the systems investigated in this 
paper. For this calculation, we assumed that the magnetic field in- 
tensity of the planet is = 14 G and of the star is S, = I G and 
that the stellar coronal material is at temperature 2 x 10^ K and has 
a base density of 10^ cm"^. 

In the second scenario investigating changes on the size of 
the planetary magnetosphere, we still considered that the stellar 
corona is steady, but allowed non-axisymmetry so that the planet 
shocks with material of different characteristics during its orbit. In 
this case, differences in the surrounding material will cause vari- 
ation on the size of the planet's magnetosphere, and therefore on 
the stand-off distance observed during transit observations. The re- 
sultant shock variation can occur both in circularised and eccentric 
systems. Because the stellar rotation period in general differs from 
the orbital period, a series of transit observations can probe differ- 
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ent stellar material. This is the case, for instance, if the star has 
an oblique magnetosphere and the planetary orbit takes the planet 
through regions of confined and expanding stellar material. 

In the third scenario investigated, we considered time- 
dependent intrinsic variations of the stellar magnetism, such as due 
to CMEs and stellar magnetic cycles. Although the impact of a 
CME on a planet increases the local density surrounding the planet, 
its effect is transient and may not be captured in transit observa- 
tions, except for the case of planets orbiting young, magnetically 
active stars for which the magnitude and frequency of flares and 
coronal mass ejections may be much higher than for the present- 
day Sun. 

Our results show that, only in the case of a circularised planet 
orbiting a perfectly aligned star (axisymmetric corona where the 
stellar rotation axis is parallel to orbital axis), a bow shock sur- 
rounding a planet should be steady. Therefore, for systems where 
a shock is detectable through near-UV transit light curve observa- 
tions, shock variations should be a common occurrence, seen in 
time-variability of near-UV transits. 
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